Spin and orbital excitation spectrum in the Kugel-Khomskii model 
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We discuss spin and orbital ordering in the twofold orbital degenerate superexchange model in 
three dimensions relevant to perovskite transition metal oxides. We focus on the particular point 
on the classical phase diagram where orbital degeneracy is lifted by quantum effects exclusively. 
Dispersion and damping of the spin and orbital excitations are calculated at this point taking into 
account their mutual interaction. Interaction corrections to the mean-field order parameters are 
found to be small. We conclude that quasi-one-dimensional Neel spin order accompanied by the 
uniform d 3z 2_ r 2-type orbital ordering is stable against quantum fluctuations. 
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It is well known that the orbital (quasi)degeneracy of 
3d-states in transition metal oxides plays an important 
role in their magnetic and lattice properties. An orbital 
ordering driven by exchange interactions and/or by Jahn- 
Teller effect occurs at low temperature resulting in a rich 
variety of magnetic structures (for a review see, e.g. [jjj). 
On the other hand a little is known, however, on dy- 
namical aspects of the coupling between spin and orbital 
degrees of freedom in these systems: a) What is the spec- 
trum of low-energy orbital excitations, b) How orbital ex- 
citations are coupled to the spin sector, c) How this cou- 
pling affects magnetic order parameter and spin waves? 
In the present paper we address these questions by con- 
sidering the superexchange model with twofold orbital 
degeneracy, which corresponds to the d 9 Mott-Hubbard 
insulator on a cubic lattice. 

To be specific, we consider the following Hamiltonian 
derived by Kugel and Khomskii Q , and studied recently 
by Feiner et al. 0]: 
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In Eq.(l) we follow notations used in g: t is the hop- 
ping between e g (3z 2 — r 2 ) orbitals along the c-axis, Si is 
the spin 1/2 operator. Operators r" act in the orbital 
subspace with basic vectors ( Q ) , (°) corresponding to the 
e g (x 2 — y 2 ) ~| x) and e g (3z 2 — r 2 ) ~| z) orbital states 
respectively. The structure of r" depends on the index 
a which specifies the orientation of the bond (ij) relative 
to the cubic axes a, b, c: 



where a z and a x are Pauli matrices. 

It is rather easy to see that the classical Neel state (i.e. 
(SiSj) = —1/4 in Eq. (1)) is infinitely degenerate: or- 
bitals at each site may be rotated independently. Feiner 
et al. H have suggested that local orbital fluctuations as- 
sociated with this degeneracy strongly affect spin-sector 
when quantum fluctuations around the Neel state are in- 
cluded, and drive the system into a disordered spin-liquid 
state even in three dimensions || . Our results presented 
below do not support this exciting scenario. We have 
investigated spin and orbital orderings, and their exci- 
tations in the model defined by Eqs. (1),(2). Our main 
findings are i) \ z)-type orbital ordering favouring quasi- 
one-dimensional spin order is the most promising candi- 
date for the ground state, ii) Orbital excitations have 
a gap generated by quantum effects. This gap controls 
well the fluctuations around the mean-field solution. Hi) 
Spin-orbit coupling does indeed act to decrease the stag- 
gered moment, but this effect is not enough to destroy 
the long-range order in a cubic lattice. 

To begin with, we use the condition Y^u j\ T t — fol- 
lowing from Eq. (2), and represent Eq. (1) in a more 
transparent way: 
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The first term in (3) represents the classical Neel en- 
ergy (in units of t 2 /U) which we drop hereafter. From the 
above Hamiltonian the key feature of the Kugel-Khomskii 
model is evident: The exchange "constant" has in fact an 
internal operator structure accounting for the orbital dy- 
namics, and it's expectation value strongly depends on 
the orientation of orbitals. It follows from Eqs. (3), (4) 
that the strength of the intersite orbital coupling (hence 
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the energy gain due to the orbital ordering) is propor- 
tional to the deviation of spins from the Neel state, i.e. 

to the value of (SiSj + —). This acts to reduce the ef- 
fective dimensionality of the spin system: Orbitals are 
arranged in such a way which makes the exchange cou- 
pling strongly nonuniform thus enhancing spin fluctua- 
tions as much as possible. In low-dimensional models, 
a similar consideration suggest that the orbital ordering 
may lead to the spin-liquid state Q] . The z-type ordering 
of orbitals in the model (1) is suggested by this picture. 
Indeed, the expectation value of exchange coupling (4) 
between z orbitals is J c = 4 along the c axis, and it is 
only small in the (ah) plane: J± = 1/4. Exchange energy 
is mainly accumulated in c chains and can be approxi- 
mated as J c (SiSj + -) c + 2J ± (S l S j + ^ -0.65 per 

site (using (SiSj) c = 1/4 — ln2 for ID Q and assuming 
(SiSj}± ~ 0). On the other hand x-type ordering results 
in the easy plane magnetic structure (J a ,b = 9/4, J c = 0) 
with a much smaller energy gain ~ —0.38. 

Our strategy is to study the Hamiltonian (3) within 
the following scheme, i) We rewrite (3) in the form 
H = H sp + H or b + H int . Here the first two terms de- 
scribe spin and orbital sectors in the mean-field level: 

= + J)» ( 5 ) 

(id) 

H orb = '£0 i S j + ±)5(JV), (6) 

(id) 

where 5 A = A — (A). The crucial importance is the 
stability of the mean-field state against fluctuations gen- 
erated by dynamical coupling between spin and orbital 
excitations. This coupling is represented by 

H mt = Y J 5{J l < l)5{S l S ] ). (7) 

(id) 

ii) We assume the antiferromagnetic spin order and uni- 
form z- or a;-type ordering of orbitals. Then we em- 
ploy spin (orbital) wave representation for Si (pi) op- 
erators. Hi) We calculate spin-orbit interaction correc- 
tions to the excitation spectrum and to the order pa- 
rameters. Since latter quantities enter in coupling con- 
stants in Eqs. (5), (6), all steps have to be done in the 
self-consistent way. 

Consider z orbital order which results in a highly 
anisotropic quasi- ID magnetic structure. We discuss first 
mean- field results, which follow from Eqs. (5), (6). Spin 
and orbital wave energies are given by u>ik = J\ \J 1 — 7ft , 
and CL>2fc = Ji\f \ + 27 2fe respectively. Here J\ = (J c + 
2J±), and 

J c = (l-a!-a*+at^) c , (8) 




The orbital stiffness is controlled by J2 = — 8(k c — —kx), 

with K a = (SiSj + -) a . Momentum dependencies of 

Lj n k (index n = 1,2) are determined by the functions 

3 k 

lik = (J c cosfc z + 2Jj_7fe)/Ji, 72fc = 7fe, 

2 4 — k 

where k = k±/k c , and 7^ = -(cosk x + cosk y ). 

We calculate all expectation values within linear spin 
(orbital) wave theory, with only one exception, the 
interchain spin correlator which we approximate as 
(S i S j ) x = (Si)(Sp + (S+Sr) |. 

Self-consistent mean-field calculations show that the 
orbital pseudospin is almost saturated (the mixture of 
I a;) state is about one percent only). Coupling between 
chains J± is weak (see Table 1) but sufficient to pro- 
duce quite large magnon dispersion in the (ab) plane (see 
thin dashed lines in Figs. 1,2). Orbital excitations are 



TABLE I. Neel order parameter (S z ) and some other ex- 
pectation values (see text for notations) calculated in the self- 
consistent mean-field approximation (V = 0), and corrected 
by including fluctuation effects (V 7^ 0). 
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0.226 


0.052 


-0.417 


-0.122 


-0.609 


-0.609 




0.191 


0.072 


-0.421 


-0.103 


-0.564 


-0.690 




FIG. 1. Dispersion ujk and damping 7^ of the spin 
and orbital waves along the direction r — > M (n, 7r, 0) — > 
R(n, 7T, it) — ► r in the Brillouin zone, calculated includ- 
ing fluctuation effects (solid lines), and in the mean-field 
approximation (dashed lines). Thin (thick) lines cor- 
respond to the spin (orbital) excitation. 7^ for orbital 
waves is almost indistinguishable from the zero line. 
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gapfull, since the orbital ordering is not associated with 
the breaking of any continuous symmetry. Of the similar 
spirit mean-field picture was recently discussed by Ishi- 
hara et al. |Q in context of their spin-orbital model for 
manganites. Quantitatively, we find that the orbital gap 
is smaller than the spin-wave bandwidth. The softness 
of the orbital excitations is related to that the orbital 
degeneracy in the model (1) can be lifted only due to 
quantum effects in the spin sector. 

2 I 1 1 1 
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FIG. 2. The same as Fig. 1 but along the T - 
X(w, 0, 0) -> N(ir, 0, tt) -> T direction. 

Now, what happens when we switch on the coupling 
between spin and orbital excitations? The latter is rep- 
resented by Hi n t (7), which in terms of spin (/?/.) and 
orbital (ipk) wave excitations reads as 

H int = -Vj2 + (9) 

k p 

\g k APiP-p + P-k0p)}(<p* + <pt q ), 

where the lowest order (three magnon) terms are only 
kept. Here V — Vo/2, and q = k — p. Matrix elements 
are 

fk,p = (U2q + V2q)[f] q Mk,p + {r)k + Vp)^k,p], (10) 

9k, P = (u2q + v 2 q)[r)qNk,p + {rjk + r)p) M k,p}, 
M k , p = (uikuip + vikVip), Nk, P = (uifcWip + V\kU\ p ), 

with rjk = (cosfcj; — cosfc a )/2. The Bogoliubov trans- 
formation coefficients in the spin subspace are given by 
uik ={(s + 1)/2} 1 /2 j Vlk = _{( s _ l)/2} 1 /2 sgn7lfcj and 

* = ( 1 -7i 2 fc ) _1/2 - The factor {u 2k + v 2k ) = (l + 2 72fc )- 1 / 4 
in (10) is due to Bogoliubov transformation in the orbital 
sector. 

Physically, the interaction (9) accounts for the process 
when spin exchange is accompanied by the simultane- 
ous orbital transition | z) «-*| x), thus enhancing the 



x orbital component in the ground state. Spin-orbital 
coupling leads to the conventional 2x2 matrix bosonic 
Green's function in both subspaces, with a diagonal (G) 
and nondiagonal (F) components given by 

G w , fc = [(iw - A u ^) + (oj k + S u , k )]/Det, (11) 
F^, k = -S^i/Dct, 
Dct = (iu - A u , k f - 

(w fc + S^ k - EjW)(w fc + S^ k + ) 

Here A u ^k and S u k represent the antisymmetric and 
symmetric (with respect to the Matsubara frequency ito) 
components of the diagonal self-energy E^l respectively, 

while S^ 1 ]. is a nondiagonal element of the self-energy 
matrix. It is implied that all quantities in Eq. (11) carry 
the subspace index n as well, and n = 1(2) stands for 
spin (orbital) waves. We calculate self-energies from the 
lowest order diagrams shown in Fig. 3. 




i) ii) Hi) 

FIG. 3. Spin-orbit interaction corrections to the spin 
(i) and orbital excitations (ii) , and to the ground state 
energy (Hi). Lines (wavy lines) represent spin (orbital) 
waves. 



In spin subspace we find (at zero temperature): 

A^.k = V 2 y^(/fc,p — 9k,p)-, — ; ■ rs 2 T~ (12) 

Here e k , p = (uJxp + ^2q), q = k — p, and 

(f - 9 2 )k,p = [Vq - (Vk + Vp) 2 ]x2q, (13) 
(/ ± g)k,p = (Vk +Vp± Vq) 2 {xikXl p ) ±1 X 2 q, 

where x lk = [(1 - 7u)/(l + 7ifc)] 1/2 , x 2k = (1 + 
2 72fc )- 1/2 . 

In the orbital sector one finds that A^^ = 0, and 

S u ,k = ^ k = 2y2 g2 +v-P ( w + J)2_ ^2 ( 14 ) 

where e k:P — (u>t p + t*)x,k+p), and ff 2 +PiP can be found 
from Eq. (13). We recall that the "bare" energies oj nk in 
Eqs. (11)-(14) are also affected by the interaction, due to 
the renormalization of parameters J n and 7„/-. 

Results of self-consistent calculations by including 
interaction corrections are presented in Table 1 and 
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Figs. 1,2. Dynamical spin-orbit coupling results in the 
following: i) It enhances quantum fluctuations in both 
subspaces thus reducing the staggered moment, and in- 
creasing the weight of the x orbital component (which 
is about 6 percent). The latter effect is also reflected 
in a larger value of the ratio J±/J c ■ ii) Spin and or- 
bital excitations are both softened, which is more pro- 
nounced in the k y = plane (and in equivalent ones), see 
Fig 2. The orbital gap still remains well defined. A spin- 
Peierls like instability is absent, because of the vanishing 
matrix elements in Eq. (9) for momenta along z (note 
77(0, 0, q z ) — 0), and because of the finite interchain cou- 
pling. Hi) Spin waves get a finite damping. Orbital waves 
are almost undamped since the density of spin states in- 
side the orbital gap is small, iv) Joint spin-orbital fluc- 
tuations significantly lower the ground state energy (see 
Table 1). The latter is given by Eq — E m f -\- (Hi n t), 
where the interaction correction to the mean- field result, 
calculated from the last diagram in Fig. 3, is 

(H m t) = -V 2 ^ 9k,p/(vik + Uip + W2,fc-ja)- (15) 

k,p 

The exchange energy gain E = —0.69 per site is found, 
which is close to our above estimation from physical con- 
siderations. Summarizing, interaction effects do not qual- 
itatively change predictions of the self-consistent mean- 
field theory, which seems to work quite reasonably. This 
is an important observation giving some credit to the 
mean-field ansatz in studying more complicated spin- 
orbital models. Of course, the latter fails when the orbital 
gap is softened close to the phase boundaries between dif- 
ferent orbitally ordered states, and the dynamical spin- 
orbit coupling becomes of crucial importance. 

Considering x-type ordering, we found it to be unsta- 
ble against fluctuations. It turns out that orbital exci- 
tations around this mean-field state are gapless at the T 
point, LL>2q ~ q- In addition, the spin-orbit interaction 
vertex remains finite at q = 0, since the orbital pseu- 
dospin is not conserved quantity, and orbital waves can 
not be considered as Goldstone modes. All these lead to 
the divergencies in perturbation theory indicating that 
an x-type ordered state is not an appropriate one, as we 
already mentioned above. This result is consistent with 

In summary, we have studied the spin-orbital coupling 
problem in the specific model, where this coupling is par- 
ticularly important because of infinite degeneracy of the 
classical Neel state in this model. The problem of the or- 
bital frustration pointed out in is actually removed by 
reducing the effective dimensionality of the spin system. 
Quantum spin fluctuations then generate an orbital ex- 
citation gap through the spin-orbit coupling mechanism. 
Orbital degeneracy in the model (1) should manifest itself 
in a strong reduction of the Neel temperature, by favour- 
ing soft quasi- ID spin structure. This is consistent with 



a basic idea of Feiner at al. [|2j that orbital degeneracy, 
in general, acts to enhance quantum spin fluctuations. 
In contrast to (2| we find, however, that this effect is 
not strong enough to destroy the Neel order. Melting 
of the long-range magnetic order by orbital fluctuations 
suggested in Q does not occur in a cubic perovskite sys- 
tems, by simple reason: A certain (model dependent) 
orbital ordering always results in the three dimensional 
(albeit very anisotropic) network of exchange interactions 
among spins. Three dimensionality of the spin sector and 
existence of the orbital gap are important factors stabi- 
lizing the Neel order. We believe that the orbital gap is 
a robust property of Mott-Hubbard insulators, which is 
related to the fact that the underlying symmetry in or- 
bital subspace is only discrete one. In a metallic state, 
doped holes can drastically change the situation, by in- 
ducing low-energy orbital fluctuations || . A study of the 
orbital melting in the Kugel-Khomskii model, driven by 
hole doping, deserves further work. 
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